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It is well known that quantum mechanics is incompatible with local realistic the- 
ories. Svetlichny showed, through the development of a Bell-like inequality, that 
quantum mechanics is also incompatible with a restricted class of nonlocal realistic 
theories for three particles where any two-body nonlocal correlations are allowed 
|31j . In the present work, we experimentally generate three-photon GHZ states to 
test Svetlichny's inequality. Our states are fully characterized by quantum state to- 
mography using an overcomplete set of measurements and have a fidelity of (84ibl)% 
with the target state. We measure a convincing, 3. 60", violation of Svetlichny's in- 
equality and rule out this class of restricted nonlocal realistic models. 

INTRODUCTION 

Quantum mechanics cannot be described by local hidden-variable (LHV) theories. This 
is the conclusion of Bell's seminal work, in which he derived a strict limit to the strength 
of correlations achievable by all LHV models that is violated by quantum predictions [Tj. 
Bell's original inequality did not allow for imperfections and thus it was not accessible to 
experimental tests. Clauser, Horne, Shimony and Holt (CHSH) addressed this issue and 
developed the CHSH inequality [2], which allowed for tests in actual experiments. Since 
then a growing number of experiments have been reported (for examples, see [31 IH El El 
[71 El El [ini ini 1121 ini 111 lig lini IITIIIH]), and the overwhelming experimental evidence 
from these tests is in favour of quantum mechanics, ruling out LHV theories. It should 
be noted that, while no loophole-free Bell test has been performed, the most significant 
potential loopholes, relating to detection efficiency and space-like separation of the choices 
of measurements settings, have both separately been closed [HI 1121 [13] • 

Both Bell's inequality and the CHSH inequality were formulated for testing the corre- 
lations between just two particles. For more than two particles, Greenberger, Horne and 
Zeilinger (GHZ) showed [I9] that a contradiction between LHV theories and quantum me- 
chanics can be seen directly in perfect correlations, as opposed to statistically in imperfect 
ones. Soon thereafter Bell-type inequalities for more than two particles were developed 
[20l[2Il[221[23l[21l[25l[26l[2Zl[2H]. Quantum predictions can violate such inequalities by an 
amount increasing exponentially with the particle number [20l |2I1 [231 [211 [2S1 l27j - 

All of the aforementioned inequalities are based on the assumption that local realism 
applies to each individual particle. Two-particle inequalities have been developed which are 
in conflict with quantum mechanics although they allow restricted, but physically motivated, 
nonlocal correlations [21] . These inequalities have recently been violated experimentally [30] ■ 

Svetlichny showed that even if one allows unrestricted nonlocal correlations between any 
two of the constituent particles in a three-particle setting one can still flnd inequalities 
violated by quantum mechanical predictions [31] • The correlations allowed by Svetlichny's 
model are strong enough to maximally violate three-partite inequalities, such as Mermin's 
[32] . which assume local realism for all particles involved. A violation of such inequalities 
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therefore can only rule out LHV theories, while a violation of Svetlichny's inequality directly 
rules out a whole class of nonlocal hidden- variable theories [321 [33l |3ll [351 Svetlichny's 
work has since been generahzed to the case of N particles |3Z1 EHl EH] • 

Experimental tests have been performed confirming the violation of the Mermin inequality 
|4Uj . the Mermin- Ardehali-Belinskii-Klyshko (MABK) inequality [H], and the cluster state 
inequality developed by Scarani et. al. [52]. For an even number of particles only, a 
sufficiently large violation of the MABK inequality also rules out partially non-local hidden- 
variable models [38]. The violation of the MABK inequality in ^Ij thus confirmed genuine 
four-particle entanglement and non-locality. The original Svetlichny inequality, however, 
remains untested. 

Here, we begin with a brief theoretical description of Svetlichny's inequality. We then 
experimentally produce high-fidelity three-photon GHZ states and characterize them via 
quantum state tomography [13]. Using these states, modulo standard loopholes [HI [121 [13], 
we experimentally demonstrate a convincing violation of Svetlichny's inequality. 

THEORY 

The two assumptions resulting in Bell-type inequalities are locality and realism as they 
were introduced by Einstein, Podolsky and Rosen [H]. We first review a straightforward 
method to derive the CHSH inequality from these two assumptions following an argument 
described by Peres ^Ij. Pairs of particles are distributed to two distant parties, A and 
B. Party A (B) can choose between two measurement settings a and a' (b and b'). For 
each measurement setting, two outcomes, -|-1 or —1, are possible. Realism assumes that the 
measurement outcomes are predetermined by some properties of the system investigated. 
These properties are known as hidden variables because they are not necessarily accessible to 
observation. The additional assumption of locality requires that the measurement outcomes 
on side A are independent of the measurement setting on side B, and vice versa. Thus for 
any given pair of particles the measurement outcomes have predetermined values a = ±1 
and a' = ±1 on side A and b = ±1 and 6' = ±1 on side B. These values identically satisfy 
the relations: 

^2 = a{b + b') + a'{b-b') = ±2 

S'^ = a'ib' + b) + a{b' -b) = ±2. (1) 

If, for example, 5*2 is averaged over many trials, the absolute value must be smaller than 2, 
which results in the CHSH inequality [2]: 

|E(a, b) + E(a, b') + E(a', b) - E(a', b') | < 2, (2) 

where the correlation, -E'(a, b), is the ensemble average {ab) over the product of measurement 
outcomes a and b for measurement settings a and b, respectively. 

This argument can be extended to three particles [38] • We will denote the particles as 
well as the measurement outcomes as a, b and c, the measurement settings as a, b, c. The 
outcome of each measurement can be +1 or —1. If we assume local realism for each of the 
three particles, then for a given set of three particles the measurement outcomes a, b and c 
as well as their primed counterparts will have predetermined values ±1. Using ([T]) we find 
that the following identity must hold: 



5*3 = 5*2(0 + c') + 5*2(0 — c) = 2 {a be + ab'c + abc — a'b'c) = ±4. 



(3) 
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Dividing this expression by two, and averaging over many trials yields Mermin's inequality 
[20] for three particles: 

|E(a', b, c) + E(a, b', c) + E(a, b, c') - E(a', b', c') | < 2, (4) 

where E{a,h,c) = (abc). 

Now assume that we allow arbitrary (nonlocal) correlations between just two of the 
particles, say a and b, while we still assume local realism with respect to the third particle, 
c. In this case we cannot factorize 5*2 as we did in ([T]) because outcomes for particle a might 
nonlocally depend on the outcomes and/or measurement settings for particle b. However, 
we can still write 

^2 = (ab) + {ab') + (a'b) - {a'b') 

S'^ = [a'b') + (a'b) + {ab') - (ab), (5) 

where the parentheses are meant as a reminder that these quantities should be regarded as 
separate and independent quantities. Each of these quantities as well as c and c' must take 
predetermined values ±1 because we assume local realism with respect to the third particle. 
This model is strong enough to violate, and reach the algebraic maximum of Mermin's 
inequality (since 5*2 can be ±4). Thus no experimental violation of Mermin's inequality can 
rule out this restricted nonlocal hidden- variable model. 

With this in mind let us slightly modify our argument to derive Svetlichny's inequality. 
Because S2 and S'2 are functions of the same four quantities, they are not independent. For 
example, whenever one of the two quantities reaches its algebraic maximum ±4, the other 
one will be 0. As a result the following identity holds: 



S2C — S2C = (ab) c + (ab) c + {ab') c — {ab') c + {a'b) c — {a'b) c — {a'b') c — {a'b') c' 

= ±4, ±2,0. (6) 

Averaging over many trials yields the Svetlichny inequality: 

= |E(a,b,c) + E(a,b,c') + ^(a,b',c) -E(a,b',c') + 

E{a', b, c) - E{a', b, c') - E(a', b', c) - E(a', b', c') | < 4, (7) 

where we refer to as the Svetlichny parameter. It is remarkable that, although we started 
out by allowing nonlocal correlations between particles a and b while c is local, one gets 
an expression identical to ([T]) if 6 and c are nonlocally correlated while a is local, or if a 
and c are nonlocally correlated while b is local. Every hidden-variable model that allows for 
nonlocal correlations between any two particles but not between all three can be seen as a 
probabilistic combination of models where the partition of the particles between nonlocal 
and local is made one or the other way. All of these models fulfill the Svetlichny inequality 

m- 

It was shown by Svetlichny that his inequality can be violated by quantum predictions, 
and that the maximum violation can be achieved with GHZ states. Assume we have a 
polarization-entangled GHZ state \ip) = {\HHV) + \VVH)), and let our measurement 
settings all be in the xy-plane of the Bloch sphere, i.e. we can write the corresponding states 
we project on as 4= {\H) + e^^lV)). For example, the measurement settings a and a' for 



4 



particle a correspond to projective measurements on the states |v4(±)) = ± e^'^''\V)) 

and |A'(±)) = ie^'^^lV^)), respectively. Here, the ± corresponds to the state the 

particle is projected on if the outcome of the measurement is ±1. For particles b and c we 
choose an analogous notation. Then the quantum prediction for the left-hand side of ([T]) is 

I cos (0a + 4>b- <Pc) + COS (0a + 06 " 0c) + COS (0a + 0fe - 0c) " 
COS (0a + 06 - 0c) + COS (0'„ + 0b - 0c) - COS (0'„ + 06 - 0c) - 

COS (0'„ + 0';, - 0c) - COS {(j)'^ + 0; - 0'J I . (8) 

With a suitable choice of angles, such as: 

0a = Y,0'a = ^,06=-,0; = O,0c = O,0^=-, (9) 

this results in = 4a/2 ^ 4, which is the maximum violation of Svetlichny's inequality 
achievable with quantum mechanics [34jj- 

Since any hidden-variable model describing a three-particle state where only two particles 
are nonlocally correlated has to fulfill the Svetlichny inequality, its violation explicitly rules 
out this type of nonlocal hidden-variable theory. 



EXPERIMENT 

Our experiment uses a pulsed titanium: sapphire laser (rep. rate 80MHz, 2.5 W avg. 
power, 790 nm center wavelength, 9nm FWHM bandwidth). We frequency double the near- 
infrared beam, producing 700 mW average power near 395 nm with a FWHM bandwidth 
of 1.8 nm. This upconverted beam is focused on a pair of orthogonally-oriented /3-Barium- 
Borate (BBO) nonlinear crystals [46] cut for type-I noncollinear degenerate spontaneous 
parametric down-conversion (SPDC) with an external half opening angle of 3°. The pump 
polarization is set to \D) = -^{\H) + \V)) such that each pump photon can produce a 
photon pair either in the first or the second BBO crystal. To compensate for temporal 
distinguishablity between the pairs created in the first and the second crystal the pump 
passes through a 1mm a-BBO crystal, a 2mm quartz crystal, and a 0.5mm quartz crystal, 
all cut for maximum birefringence. To compensate the 75//m spatial walk-off between the 
horizontally and vertically polarized SPDC photons observed in one of the output modes, 
we insert a 0.75mm thick BiBO crystal cut at 6 = 152.6° and = 0°. For these cut 
angles the crystal compensates the transverse walk-off without introducing additional time 
walk-off. The photons are subsequently coupled into single-mode fibres. We label the two 
corresponding spatial output modes as 1 and 2, see figure 1(a). Fiber polarization controllers 
ensure that in mode 1 states in the HV basis remain unchanged while in mode 2 we flip the 
polarization, i.e., H ^ V. With this configuration we achieve a two-photon coincidence rate 
of 43 kHz and single rates of about 240 kHz and 270 kHz for modes 1 and 2, respectively. 
The measured contrast of the pairs is 75 : 1 in the H/V basis and 61 : 1 in the 45°/ — 45° 
basis when the source is adjusted to produce |0"'") = i\HH) + |V^V)) states. 

Following the approach in [17] we use the double-pair emission of the SPDC source to 
produce 3-photon GHZ correlations in the interferometer shown in figure 1(b). A four- 
fold coincidence detection in the four outputs of the interferometer indicates the successful 
generation of the GHZ state. To lowest significant order, a four-fold coincidence can only 
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FIG. 1: Experimental setup, (a) Schematic of our type-I SPDC source. A 45° polarized, pulsed 
UV beam pumps a pair of orthogonally oriented BBO crystals cut for type-I phase matching. 
Temporal walk-off between the pairs created in the first and in the second crystal is compensated 
by a 1 mm thick q-BBO and two quartz crystals (one 2 mm and one 0.5 mm thick) before the 
SPDC crystals. Spatial walk-off, which occurred in mode 1, was compensated by a BiBO crystal. 
The phase between horizontal and vertical photons was adjusted by tilting a A/4 plate in mode 
2. All photons pass through 3nm FWHM bandwidth filters around 790 nm and are coupled into 
single-mode fibres corresponding to the spatial modes 1 and 2. (b) The interferometer used to 
project on a three-photon GHZ state. Using fibre polarization controllers, the polarization in 
mode 2 is rotated such that we map H V and V ^ H, while in mode 1, H and V are preserved. 
Inside the interferometer the four photons from a double-pair emission can be split up into four 
separate spatial modes and result in a four- fold coincidence event between the detectors T, Da, 
Db, and D^. In this case the three photons in the modes a, b, and c will be projected on the 
three-photon GHZ state ^ {\HaHbVc) + \ VaVbHc)) given that the photons detected by Da and Db 
are indistinguishable. 

occur if two photons enter the interferometer via the spatial mode 1, and their two partner 
photons enter it via mode 2. The two photons in input mode 1 impinge on a polarizing 
beam splitter (PBS). In order for a four-fold coincidence event to occur one of them must be 
H polarized and the other V polarized. The V photon is reflected at the PBS. Its detection 
by detector T serves as a trigger event. The H photon passes through the PBS in mode 1, 
and a A/2 plate oriented at 22° rotates the polarization from \H) to \D). The two photons 
in mode 2 are split at the 50 : 50 beam splitter (BS) with probability 1/2. Only in this case 
a four-fold coincidence can occur. The transmitted photon and the \D) polarized photon in 
mode 1 are overlapped on a PBS. A coincidence detection event in the two output modes 
of the PBS can only occur if both photons are transmitted or both are reflected. If these 
two possibilities are indistinguishable, Hong-Ou-Mandel interference will occur [58]. In the 
reflected mode of the BS we compensate for a phase shift due to birefringence in the BS by 
tilting a A/4 plate. 

A four-fold coincidence detection in the interferometer outputs can only occur if the 
trigger photon is V polarized, and if the other three photons in the modes a, b and c (see 
figure 1(b)) are either \HaHbVc) or \VaVbHc). By tilting a A/4 plate in output mode a we 
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FIG. 2: Reconstructed three- photon density matrix, a) Real part and b) imaginary part of the 
density matrix. The state was reconstructed from a tomographically overcomplete set of 216 
measurements. For each measurement the four- fold coincidence counts were integrated over 27min, 
see table ll The fidelity of the reconstructed density matrix with the GHZ state -^{\HHV) + 
\VVH)) is (84±1)%. 

adjust the relative phase between these contributions such that a four- fold coincidence event 
signals a GHZ state of the form ^ {\HaHbVc) + \VaVbHc)). 

Given this state, quantum mechanics predicts a maximum violation of Svetlichny's in- 
equality if we choose measurements of the form \H) + e^'^\V) with the angles given in 
Particles a, b and c are identified with photons in the interferometer output modes a, b and 
c, respectively. To test Svetlichny's inequality each of the photons has to be measured in 
two measurement bases, and for each basis there are two possible outcomes, +1 and — 1. 
For each outcome we have to set the polarization analyzer in the respective mode such that 
a photon detected after passing through the analyzer corresponds to that outcome. 

The Svetlichny parameter, (see ([T])), consists of 8 correlations, each of which can be 
constructed from 8 three-photon polarization measurements for a total of 64 measurements. 
Each polarization analyzer consists of a A/2-plate followed by a A/4-plate and then a PBS. 
The photons passing the PBSs are detected with single-mode fibre-coupled single-photon 
counting modules; the coincidence window is 10 ns. 

In order to fully characterize the state produced by our setup we perform quantum state 
tomography. Because all the measurement settings for the Svetlichny inequality lie in the 
xy-plane of the Bloch sphere these alone are not tomographically complete. Instead of 
performing an additional run apart from the measurements of the Svetlichny settings, we 
add two additional projective measurements {\H) and |V^)) for each of the particles, resulting 
in a total of 216 three-photon polarization measurements. Our set of measurements is now 
actually tomographically overcomplete, which has been shown to produce better estimates 
of quantum states |19] . 
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TABLE I: Experimentally-measured counts. Four-fold coincidences for the 216 measurements 
performed for quantum state tomography. A subset of 64 measured counts were used to test 
Svetlichny's inequality; these counts are shown in boldface type. We cycled through all of the 
measurements 27 times, counting for 60s for each measurement. The four-fold coincidences given 
are the result of integrating over all of these cycles. 

RESULTS 

In our setup the maximum four-fold coincidence rates of 7 x 10~^ Hz and 8 x 10~^ Hz 
were achieved for the correlations \HaHbVc) and \VaVbHc), respectively. To get statistically 
significant counts we integrated over 27min per measurement. In order to reduce the nega- 
tive effects of misalignment of the setup over time we realized this integration by counting 
for 60 s for each of the 216 measurements, then repeating the full cycle of measurements 27 
times. The resulting counts are given in table [T| 

We applied the maximum likelihood technique pOj to reconstruct the density matrix of 
our state. Its real and the imaginary part are shown in figure 2. The fidelity F = {ip\p\ip) 
of the density matrix with the GHZ state ^ {\HHV) + \VVH)) is 0.84 ± 0.01. The errors 
of quantities derived from the reconstructed density matrix were calculated via a Monte 
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FIG. 3: Measured correlations and Svetlichny parameter, a) Measured correlations for the eight 
combinations of measurement settings for the three particles. Each correlation is constructed 
from 8 four-fold coincidence measurements integrated over 27min. The count rates for each of 
these measurements are given in table |l| b) These correlations yield a Svetlichny parameter of 
4.51 lb 0.14, which clearly violates the bound (dashed line) of 4 of the Svetlichny inequality. The 
quantum mechanical limit is 4\/2. Even higher values (pattern-filled region) can be reached by 
allowing arbitrarily strong nonlocal correlations. 



Carlo simulation, where we used each of the measured counts as the mean of a Poissonian 
distribution. According to these distributions we generated random counts and ran the 
maximum likelihood algorithm. This procedure was repeated 400 times, and we report the 
standard deviation and mean for quantities derived from these reconstructed states. 

The 64 measurements that quantum mechanics predicts to violate Svetlichny inequality 
are among the 216 measured. After integrating over all 27 cycles we get a Svetlichny pa- 
rameter of Sy = 4.51 ± 0.14; the eight measured correlations are shown in figure 3. This 
value violates the Svetlichny inequality by 3.6 standard deviations. It is, however, in good 
agreement with the value, iS^^^ = 4.48 ± 0.11, predicted by quantum mechanics given the 
reconstructed density matrix. 



CONCLUSION 

We used the double-pair emission from a pulsed type-I SPDC source and projected the 
photons onto a GHZ state using a linear optical interferometer. We fully characterized 
the generated state and reconstructed the density matrix applying the maximum likelihood 
technique [50] using an overcomplete set of measurements. From the reconstructed density 
matrix, we found that our state matched the target GHZ state with a fidelity of (84 ± 1)%. 

We experimentally demonstrated the violation of the original Svetlichny inequality for 
a three-particle GHZ state with a value of 4.51 ± 0.14, which is greater than 4 by 3.6 
standard deviations. This value is in good agreement with that predicted by quantum 
mechanics from our reconstructed density matrix, 4.48 ± 0.11. By violating Svetlichny's 
long-standing inequality, we have shown that the correlations exhibited by three particles 
cannot be described by hidden-variable theories with at most two-particle nonlocality. 
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